Production of a sterile species: quantum kinetics. 



D. Boyanovskyi-B and C. M. Ho^'Q 

^Department of Physics and Astronomy, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA 

(Dated: February 1, 2008) 

Production of a sterile species is studied within an effective model of active-sterile neutrino mixing 
in a medium in thermal equilibrium. The quantum kinetic equations for the distribution functions 
and coherences are obtained from two independent methods: the effective action and the quantum 
master equation. The decoherence time scale for active-sterile oscillations is r^ec = '2/'^aa, but the 
evolution of the distribution functions is determined by the two different time scales associated with 
the damping rates of the quasiparticle modes in the medium: Fi = Faa cos'^ Om \ T2 = F^a sin'^ 9m, 
where Fan is the interaction rate of the active species in absence of mixing and 6m the mixing angle 
in the medium. These two time scales are widely different away from MSW resonances and preclude 
the kinetic description of active-sterile production in terms of a simple rate equation. We give the 
complete set of quantum kinetic equations for the active and sterile populations and coherences 
and discuss in detail the various approximations. A generalization of the active-sterile transition 
probability in a medium is provided via the quantum master equation. We derive explicitly the 
usual quantum kinetic equations in terms of the "polarization vector" and show their equivalence 
to those obtained from the quantum master equation and effective action. 
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I. INTRODUCTION 



Sterile neutrinos, namely weak interaction sing lets, are ac quir ing renewed attention as potential candidates for cold 
or warm dark matter P, H a li, H, 0, H, S Efl IH [13, [Slli, and may also be relevant in stellar collapse [l^ ITgI . 
primordial nucleosvnthesis|17l Il8l| . and as potential explanation of the anomalous velocity distributions of pulsars [l9|, 
[20, 21]. Although sterile neutrinos are ubiquitous in extensions of the standard model[22, 23, 24, 25], the MiniBooNE 
collaboration[2g has recently reported results in contradiction with those from LSNDj27l. |28| that suggested a sterile 
neutrino with Am^ ~ 1 eV^ scale. Although the MiniBooNE results hint at an excess of events below 475 MeV the 
analysis distinctly excludes two neutrino appearance-only from v^^ — > oscillations with a mass scale Am^ ~ 1 eV^, 
perhaps ruling out a light sterile neutrino. However, a recent analysis[2^ suggests that while (3 -I- 1) schemes are 
strongly disfavoured, (3 -I- 2) neutrino schemes provide a good fit to both the LSND and MiniBooNE data, including 
the low energy events, because of the possibility of CP violation in these schemes, although significant tension remains. 

However, sterile neutrinos as dark matter candidates would require masses in the keV range[l|, d, H, [3, [^, 0, 0, S 
[Til . [T3 ]. hence the MiniBooNE result does not constrain a heavier variety of sterile neutrinos. The radiative decay 
of keV neutrinos would contribute to the X-ray background [1, [10] ■ Analysis from the X-ray background in clusters 
provide constraints on the masses and mixing angles of sterile neutrinos [0, [31I [3^ . [33j . and recently it has been 
suggested that precision laboratory experiments on (3 decay in tritium may be sensitive to ^ keV neutrinos [3^. Being 
weak interaction singlets, sterile neutrinos can only be produced via their mixing with an active species, hence any 
assessment of the possibility of sterile neutrinos as dark matter candidates or their role in supernovae must begin 
with understanding their production mechanism. Pioneering work on the description of neutrino oscillations and 
decoherence in a medium was cast in terms of kinetic equations for a fiavor "matrix of densities" 35] or in terms 
of 2 X 2 Bloch-type equations for fiavor quantum mechanical stat es l36l . [37j . A general field theoretical approach 
to neutrino mixing and kinetics was presented in [38l . [3^ (see also |25|]). however, while such approach in principle 
yields the time evolution of the distribution functions, sterile neutrino pro duction in the early Universe is mostly 
studied in terms of simple phenomcnological rate equationsP, [3, [13, [4l[. l42i. i43j . An early approach ^40] relied on a 
Wigner-Weisskopf effective Hamiltonian for the quantum mechanical states in the medium, while numerical studies of 
sterile neutrinos as possible dark matter candidates P. [43| rely on an approximate approach which inputs an effective 
production rate in terms of a time averaged transition probability 41, 42]. More recently the sterile production rate 
near an MSW resonance including hadronic contributions has been studied in ref . (5] . 
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The rich and complex dynamics of oscillations, decoherence and damping is of fundamental and phenomenological 
importance not only in neutrino cosmology but also in the dynamics of meson mixing and CP violation^45l l46j . In 
ref.fi^] it was argued that the spinor nature of neutrinos is not relevant to describe the dynamics of mixing and 
oscillations at high energy which can then be studied within a (simpler) quantum field theory of meson degrees of 
freedom. 

Recently we reported on a study (isjl of mixing, decoherence and relaxation in a theory of mesons which provides 
an accurate description of similar phenomena for mixed neutrinos. This effective theory incorporates interactions 
that model the medium effects associated with charge and neutral currents for neutrinos and yield a robust picture 
of the non-equilibrium dynamics of mixing, decoherence and equilibration which is remarkably general. The fermion 
nature of the distributions and Pauli blocking effects can be simply accounted for in the final result [1^. This study 
implemented quantum field theory methods to obtain the non-equilibrium effective action for the "neutrino" degrees 
of freedom. The main ingredient in the time evolution is the full propagator for the "neutrino" degrees of freedom in 
the medium. The complex poles of the propagator yield the dispersion relation and damping rates of quasiparticle 
modes in the medium. The dispersion relations are found to be the usual ones for neutrinos in a medium with the 
index of refraction correction from forward scattering. For the case of two flavors, there are two damping rates which 
are widely different away from MSW resonances. The results of this study motivated [49l a deeper scrutiny of the rate 
equation which is often used to study sterile neutrino production in the early Universe[J, |4^ . |43| . 

One of the observations in [4^ is that the emergence of two widely different damping time scales precludes a reliable 
kinetic description in terms of a time averaged transition probability suggesting that a simple rate equation to describe 
sterile neutrino production in the early Universe far away from MSW resonances may not be reliable. 

Motivation and goals: The broad potential relevance of sterile neutrinos as warm dark matter candidates in 
cosmology and their impact in the late stages of stellar collapse warrant a deeper scrutiny of the quantum kinetics of 
production of the sterile species. Our goal is to provide a quantum field theory study of the non-equilibrium dynamics 
of mixing, decoherence and damping and to obtain the quantum kinetic equations that determine the production of 
a sterile species. We make progress towards this goal within a meson model with one active and one sterile degrees 
of freedom coupled to a bath of mesons in equilibrium discussed in ref. ^4§\ . As demonstrated by the results of ref . [4^ 
this (simpler) theory provides a remarkable effective description of propagation, mixing, decoherence and damping 
of neutrinos in a medium. While ref. [4^ studied the approach to equilibrium focusing on the one body density 
matrix and single quasiparticle dynamics, in this article we obtain the non-equilibrium effective action, the quantum 
master equation and the complete set of quantum kinetic equations for the distribution functions and coherences. We 
also establish a generalization of the active-sterile transition probability based on the quantum master equation. In 
distinction with a recent quantum field theory treatment [i^l we seek to understand the quantum kinetics of production 
not only near MSW resonances, at which both time scales concidejH, but far away from the resonance region 
where the damping time scales are widely separated (isl l49j . 

Similarities and differences: The scalar field model that we study has many similarities with the neutrino case 
but also important differences. Similarities: as demonstrated in our previous studv[48t a) the scalar model describes 
"flavor mixing" in a similar manner as in the case of neutrinos, where mixing arises from off-diagonal mass matrix 
elements, b) a medium induced "matter potential" which arises from the forward-scattering contribution to the (self) 
energy, c) the dispersion relation for the propagating modes is identical to those of neutrinos in a medium [49l. [sT 
the effective mixing angles in the medium have a functional form identical to those for neutrinos in a medium[25i 
e) the form of the transition probability for ensemble averages in the medium is identical to that for the active-sterile 
neutrino transition probability [1^, f) the relationship between the damping rates of the propagating modes and the 
active collision rate is identical to the neutrino case[49l| and g) as shown in detail in section (jVI[) the kinetic equations 
obtained are identical to those in terms of the polarization vector often quoted in the neutrino literature (see section 
(jVip V Differences: There are obvious differences with the neutrino case that should not be overlooked: a) spinor 
and chirality structure: although this is a clear difference, it is important to highlight that neither the quantum 
mechanical description of neutrino mixing nor the phenomenological description of neutrino kinetics account for 
either spinorial structure or chirality. b) Fermionic vs. Bosonic degrees of freedom, the most obvious difference 
is in the distribution functions, however the results obtained in sections piUIVip for the kinetic description allow a 
straightforward replacement of the distribution functions for the Fcrmi-Dirac expressions thus automatically including 
Pauli blocking, c) An important difference is the matter potential, in the scalar model this is given by the one-loop 
Hartree self-energy, which is manifestly positive, whereas in the case of neutrinos the matter potential features a CP- 
odd and a CP-even contribution [sO] and it can feature either sign. The existence of an MSW resonance hinges on the 
sign of the self energy, in particular on the CP-odd component. However, this important difference notwithstanding, 
our study does not rely on or require a specific form of the matter potential, only the fact that the matter potential is 
diagonal in the flavor basis and in the case under consideration only the active-active matrix element is non-vanishing. 
Whether or not there is an MSW resonance depends on the specific form of the matter potential and in the case of 
neutrinos, on the CP-odd (lepton and baryon asymmetry) component of the background. Our study addresses all 
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possible cases quite generally without the need to specify the sign (or any other quantitative aspect) of the matter 
potential. 

Summary of results: 

• i: We obtain the quantum kinetic equations for production by two different but complementary methods: 

a) : the non-equilibrium effective action obtained by integrating out the "bath degrees of freedom" . This 
method provides a non-perturbative Dyson-like resummation of the self-energy radiative corrections and leads 
to the full propagators in the medium. This method makes explicit that the "neutrino" propagator in the 
medium along with the generalized fluctuation-dissipation relation of the bath in equilibrium are the essential 
ingredients for the kinetic equations and allows to identify the various approximations. It unambiguously reveals 
the emergence of two relaxation time scales associated with the damping rates of the propagating modes in the 
medium Fi = Taa cos^ 6m ; r2 = Taa sin^ 6m where Taa is the interaction rate of the active species and 9m the 
mixing angle in the medium, confirming the results of references ^48v. .49,] . These time scales determine the kinetic 
evolution of the distribution functions and coherences. 

b) : the quantum master equation for the reduced density matrix, which is obtained by including the lowest 
order medium corrections to the dispersion relations (index of refraction) and mixing angles into the unperturbed 
Hamiltonian. This method automatically builds in the correct propagation frequencies and mixing angles in the 
medium. 

From the quantum master equation we obtain the kinetic equations for the distribution functions and coherences. 
These are identical to those obtained with the non-equilibrium effective action to leading order in perturbative 
quantities. After discussing the various approximations and their regime of validity we provide the full set 
of quantum kinetic equations for the active and sterile production as well as coherences. These are given by 
equations ()4. 41114.46]) in a form amenable to numerical implementation. We show that if the initial density 
matrix is off-diagonal in the basis of the propagating modes in the medium, the off-diagonal coherences are 
damped out in a decoherence time scale Tdec = Z/Faa. The damping of these off diagonal coherences leads to 
an equilibrium reduced density matrix diagonal in the basis of propagating modes in the medium. 

• ii: We elucidate the nature of the various approximations that lead to the final set of quantum kinetic equations 
and discuss the interplay between oscillations, decoherence and damping within the realm of validity of the 
perturbative expansion. 

• iii: We introduce a generalization of the active-sterile transition probability in the medium directly based on the 
quantum density matrix approach. The transition probability depends on both time scales 1/Fi, l/r2, and the 
oscillatory term arising from the interference of the 1, 2 modes in the medium is damped out on the decoherence 
time scale Tdec but this is not the relevant time scale for the build up of the populations or the transition 
probability far away from an MSW resonance. 

• iv: We derive the quantum kinetic equation for the "polarization vector" often used in the literature directly 
from the kinetic equations obtained from the quantum master equation under a clearly stated approximation. 
We argue that the kinetic equations obtained from the quantum master equation exhibit more clearly the time 
scales for production and decoherence and reduce to a simple set within the regime of reliability of perturbation 
theory. We discuss the shortcomings of the phenomenological rate equation often used in the literature for 
numerical studies of sterile neutrino production. 

In section (jlT]) we introduce the model, obtain the effective action, and the full propagator from which we extract 
the dispersion relations and damping rates. In section (jIII[) we define the active and sterile distribution functions and 
obtain their quantum kinetic non-equilibrium evolution from the effective action, discussing the various approxima- 
tions. In section (jIV[) we obtain the quantum Master equation for the reduced density matrix, also discussing the 
various approximations. In this section we obtain the full set of quantum kinetic equations for the populations and 
coherences and show their equivalence to the results from the effective action. In section ([V]) we study the kinetic 
evolution of the off-diagonal coherences and introduce a generalization of the active-sterile transition probability in 
a medium directly from the quantum master equation. In section (jVip we establish the equivalence between the 
kinetic equations obtained from the quantum master equation and those most often used in the literature in terms 
of a "polarization vector" , along the way identifying the components of this "polarization vector" in terms of the 
populations of the propagating states in the medium 1,2 and the coherences. While this formulation is equivalent 
to the quantum kinetic equations obtained from the master equation and effective action, we argue that the latter 
formulations yield more information, making explicit that the fundamental damping scales are the widths of the 
quasiparticle modes in the medium and allow to define the generalization of the transition probability in the medium. 
We also discuss the shortcomings of the phenomenological rate equations often invoked for numerical studies of sterile 
neutrino production. Section (|VII[) summarizes our conclusions. Two appendices elaborate on technical aspects. 
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II. THE MODEL, EFFECTIVE ACTION, AND DISTRIBUTION FUNCTIONS. 

We consider a model of mesons with two flavors a, s in interaction with a " vector boson" and a "flavor lepton" 
here denoted as W, Xa respectively, modeling charged and neutral current interactions in the standard model. This 
model has been proposed as an effective description of neutrino mixing, decoherence and damping in a medium in 
ref.fUl to which we refer reader for details. As it will become clear below, the detailed nature of the bath fields \a 
is only relevant through their equilibrium correlation functions which can be written in dispersive form. 

In terms of the field doublet 



the Lagrangian density is 



r (2.1) 



^ = 1 {d^'^^d^^ - $^m2$} + Co [W, x]+GW (t>aXa + G4>lxl (2.2) 



2 

where the mass matrix is given by 

and £o[W^iX] is the free field Lagrangian density for W,x which need not be specified. 

The mesons 0a, s play the role of the active and sterile flavor neutrinos, Xa the role of the charged lepton associated 
with the active flavor and W a charged current, for example the proton-neutron current p^^{l — gAjb)'^ or a similar 
quark current. The coupling G plays the role of Gp. The interaction between the "neutrino" doublet and the W,Xa 
fields is of the same form as that studied in ref.[2^[38, 39] for neutral and charged current interactions. 

The last term in the Lagrangian density (|2.2p allows to model the matter effective potential from forward scattering 
in the medium by replacing xi by its expectation value in the statistical ensemble, (xi)- The resulting term G(j)^{Xa) 
effectively models a matter potential from forward scattering in the medium psl]. While in the bosonic case {xl} is 
manifestly positive, in the fermionic case the effective potential from forward scattering in the medium features two 
distinct contributions [50| : a CP odd contribution which is proportional to the lepton and baryon asymmetries, and a 
CP even contribution that only depends on the temperature. However, as it will become clear below, we do not need 
to specify the precise form of the matter potential or of the bath degrees of freedom, only the fact that the matter 
potential is diagonal in the flavor basis with only entry in the a — a component, and the spectral properties of the 
correlation function of bath degrees of freedom are necessary. 

The flavor ipa.s and the mass basis fields (pi.2 are related by an orthogonal transformation $ = U{d) cp 

where the orthogonal matrix U{d) diagonalizes the mass matrix M^, namely 

u-\e)M'u{e) = (^^f ^j,^ (2.5) 

In the flavor basis the mass matrix M can be written in terms of the vacuum mixing angle 9 and the eigenvalues of 
the mass matrix as 

2 V sm2t^ cos 26' / ^ ' 



where we introduced 
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M = - {Ml + M|) ; 5M^ = M| - Ml . (2.7) 

For the situation under consideration with keV sterile neutrinos with smaU vacuum mixing angle 6 

Maa ~ Ml ; Mss ~ M2 (2.8) 

and in the vacuum 

0a ~ 01 ; 4>s 4>2- (2.9) 

We focus on the description of the dynamics of the "system fields" (pa, a — a,s. The strategy is to consider the time 
evolved full density matrix and trace over the bath degrees of freedom W. It is convenient to write the Lagrangian 
density (|2.2p as 

£[<Pa,Xa,W] = £oM + Comx] + GcPaOa + Gcgxl (2.10) 

where 

Oa=XaW. (2.11) 

and £o[' • • ] are the free Lagrangian densities for the fields (f>a,Xa,W respectively. The fields <j)a are considered as the 
"system" and the fields Xa,W are treated as a bath in thermal equilibrium at a temperature T = 1//3. We consider 
a factorized initial density matrix at a time = of the form 

p(0) = p<,(0)®pb(0) ; ps(0) = e-/'^°[X''^] (2.12) 

where Ho[x, W] is Hamiltonian for the fields Xa,W in absence of interactions with the neutrino field (j)a- 

Although this factorized form of the initial density matrix leads to initial transient dynamics, we are interested in 

the long time dynamics, in particular in the long time limit. 

The bath fields Xa , W will be "integrated out" yielding a reduced density matrix for the fields (j)a in terms of an 

effective real-time functional, known as the influence functional (s^ in the theory of quantum brownian motion. The 

reduced density matrix can be represented by a path integral in terms of the non-equilibrium effective action that 

includes the influence functional. This method has been used extensively to study quantum brownian motion[53.[53l|. 

and quantum kinetics [s^. [55| and more recently in the study of the non-equilibrium dynamics of thermalization in a 

similar model [i^. The time evolution of the initial density matrix is given by 

p{t) = e-'-^(*-*»)p(to)e'-^(*-*°) (2.13) 

Where the total Hamiltonian H is 

H = Ho[<f>]+Ho[x,W]+Hi[^,X,W] (2.14) 
Denoting all the flelds collectively as X to simplify notation, the density matrix elements in the field basis are given 

by 



{X\p{t)\X') = J DX,DX[ {X\e-'"^'-'^>^\X,) {X,\p{t^)\X'^ {X'y('''"^\X') . (2.15) 
The density matrix elements in the field basis can be expressed as a path integral by using the representations 

^^|g-.i/(t-to)|^^) = yp^+ g^/4'it/';'-£[^+] . x+{tQ)^x,-x+{t) = x. (2.16) 

Similarly 

(X,'|e^^(*-*«)|X') VX- ^-^n,dtjd^xc[x-] . x-{tQ)^X[-X-{t)^X' . (2.17) 
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Therefore the full time evolution of the density matrix can be systematically studied via the path integral 

Z = J VX+VX-e'^'^o dtjd'x{c[x^]-c[x-]} ^ ^2.18) 

with the boundary conditions discussed above. This representation allows to obtain expectation values or correlation 
functions C(X, Xi, X[, X'\ t, t') which depend on the values of the fields Xi\ X[ through the initial conditions. In order 
to obtain expectation values or correlation functions in the full time evolved density matrix, the results from the path 
integral must be averaged in the initial density matrix p{to), namely 

{C{X,X';t,t')) = J DX,DX[ {X,\p[to)\X'^ C{X, X,, X^ X';t,t') . (2.19) 

We will only study correlation functions of the "system" fields (f>on therefore we carry out the trace over the Xa 
and W degrees of freedom in the path integral (|2.18p systematically in a perturbative expansion in G. The resulting 
series is re-exponcntiated to yield the non- equilibrium effective action and the generating functional of connected 
correlation functions of the fields (pa- This procedure has been explained in detail in references [s^ [55| and more 
recently in [48l | within a model similar to the one under consideration. Following the procedure detailed in these 
references we obtain the non-equilibrium effective action up to order and quadratic in the fields neglecting 
higher order non-linearities, 



effl 



I dt I dt' 



(2.20) 



where the matter potential is 



V 



(2.21) 



with the average in the initial bath density matrix. In the bosonic model, the corresponding one-loop diagram at 
order G that yields the matter potential and effectively models forward scattering in the medium is depicted in figure 
©■ 



FIG. 1: One loop self-energy for the active species at order G, corresponding to the matter potential Vaa ~ G(x^)- 



In the fermionic theory, the matter potential in a medium at finite temperature and density has two distinct 
contributions [soj : a CP-odd term proportional to the lepton and baryon asymmetries and a CP-even term that only 
depends on the temperature. The sign of these contributions may be either positive or negative depending on which 
term dominates [soj . The presence of an MSW resonance in the medium depends crucially on the CP-odd contribution. 
In the case of sterile neutrinos with masses in the keV range, only for non-vanishing lepton asymmetry is there an 
MSW resonance. However, the only important point for the analysis that follows is that the matter potential is 
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diagonal in the flavor basis, with only entry Vaa, namely the form of the matrix given by eqn. (|2.21|) . but of course 
the matrix element Vaa itself will be different for fermions. 

The correlation functions Q{t,t') ^ {Oa{t)Oa{t')) = {W{t)W{t')) {Xa{t)Xa{t')) are also determined by averages in 
the initial equilibrium bath density matrix and their explicit form is given in reference [48| (see also appendix (jB]) ). 

Performing the trace over the bath degrees of freedom the resulting non-equilibrium effective action acquires a 
simpler form in terms of the Wigner center of mass and relative variables (isi [54. Issjl 

'^a{x,t)^^{^+{x,t) + cj)-{x,t)) ; i?,(f,i) = (0+(f,t)-0-(f,<)) ; a^a,s (2.22) 

and a corresponding Wigner transform of the initial density matrix for the fields. See ref.jJsj for details. The 
resulting form allows to cast the dynamics of the Wigner center of mass variable as a stochastic Langevin functional 
equation, where the effects of the bath enter throu gh a dissipative kernel and a stochastic noise term, whose correlations 
obey a generalized fluctuation-dissipation relation]48l. [s^ [ssj . In terms of spatial Fourier transforms the time evolution 
of the center of mass Wigner fleld \E' is given by the following Langevin (stochastic) equation (see derivations and 
details in refs . [H, [13, [H [Hi, [H ) 

+ (fc' S^P + ^Ip + Va/j)* ^ git) + f dt' ^^p{k- 1 - t')^ g(i') = j:{t) 

Jo 

= 0) = ^Ij: ; i'^j:{t = 0) = n° (2.23) 

where 'i>'^ If" - are the initial values of the field and its canonical momentum. The matter potential Y in the 

equation of motion (|2.23l) effectively models the general form of the matter potential in the fermionic case. The 
specific value and sign of Vaa is not relevant for the general arguments presented below. 

The stochastic noise g(t) is described by a Gaussian distribution function [H, [H, HI] with 

/OO 7 
:5^e-(*-*')^„Mfc;^) (2.24) 

and the angular brackets denote the averages with the Gaussian probability distribution function, determined by the 
averages over the bath degrees of freedom. The retarded self-energy kernel has the following spectral representationfZsj 

/OO 
e''^(*-*')lmE„^(fc;cj)dtj (2.25) 
-OO 

where the imaginary part in the flavor basis is 

ImE(fc; u;) = ImE„„(fc; co) J [J ^ , (2.26) 

and ImSaaik','-^) is obtained from the cut discontinuity in the one-loop diagram in figure 0. In this figure the W 
propagator should be identified with the full charged vector boson propagator in the standard model, including a 
radiative self-energy correction from a quark, lepton or hadron loop. 



FIG. 2; One loop self-energy for the active species to order G^. The cut discontinuity across the W — x lines yields the imaginary 
part lmSaa{k;uj). 



Because the bath fields are in thermal equilibrium, the noise correlation kernel ICapik^oj) in eqn. (|2.24p and the 
absorptive part of the retarded self energy lmTiap{k; lo) obey the generalized fluctuation dissipation relation [isl [s^. [ssj 

'I3uj' 



K,aj3{k\ uj) — ImEQ,^(A;; uj) coth 



(2.27) 
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The solution of the Langevin equation (I2.23P is|48l. 154 l5a| 

- GMk; t) ^I j, + G^p{k- 1) nl^^^ + / G^0{k- 1') ipj:{t - t')dt' 



(2.28) 

from which is clear that the propagator Gct/3 contains all the relevant information for the non-equilibrium dynamics. 
In the Breit-Wigner (narrow width) approximation, the matrix propagator G(fc; t) in the flavor basis is given bv[48l| 



G{k- 1) = Zl'' e 



ih{k) 2 ih[k) 



where Zi" are the residues at the quasiparticle poles and wc have introduced the matrices 



(2.29) 



Ri^HO) - ^ , -cos^sinf? ^ 10 1 

^ ^ \ - cos 6* sm 6* sm^ 6* / ^ M / ^ ^ 



(2.30) 



]R(2)(6l) = 



''f\ cos0sine \ [^-i(^?) 

cos 9 sm fc* cos^ fc* / ^ M 1 / ^ ' 



(2.31) 



Il(3)(6l) ^ sin 261 



sin 26* cos 29 
cos 26* -sin 26* 



sin2(?C/(0) ( 5 J ) U-\9). 



(2.32) 



From the results of reference |48l| to leading order in G, the mixing angle in the medium is determined from the 
relations 



cos 29 — TTft 
cos 26*™ = ^ 



sin 29„ 



sin 26* 



where 



(cos2^-^)V(sin2^)^ 



(2.33) 



(2.34) 



The expressions (I2.33P for the mixing angles in the medium in terms of the mixing angle in the vacuum and the 
matter potential is exactly of the same form as in the case of (fermionic) neutrinos in a medium [2^. [23l [2^ [2^ . An 
MSW resonance occurs whenever [H [H, H, [H 



Vaa = SM^ cos 26*. 



The propagating frequencies and widths are given bv|48l| 

ni{k) = LOiik) + Auji{k) ; ri(/c) = 



lmT,aa{k;uJi{k )) ^ 
wi(fc) 



cos 9r, 



where 



^2{k) = a;2(fc) + Ac.2(fc) ; T,{k) = I^^^<^°(fc;/;^2(fc)) ^.^-.^^^ 



wi(fc) = k +M +— — 



(2.35) 

(2.36) 
(2.37) 

(2.38) 
(2.39) 
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are the propagating frequencies (squared) in the medium including the matter potential at order G, namely the index 

of refraction arising from forward scattering, with M ■,SM'^ defined in equation (|2.7p . The second order frequency 
shifts are 



Awi(fc) 



Aw2(fc) = - 



COS^ 6jn 

27rwi(A:) 
sin^ 6>„i 



dujV 



uj — uji{k) 
lmT,aa{k\uj) 

LO - UJ2{k) 



andkS 



7(fc) = 



lrnEaa{k;u!{k)) _^ 
ujUk)-ujf{k) 



(fc) = yk^ + M 



(2.40) 
(2.41) 

(2.42) 



The relationship between the damping rates ri.2 and the imaginary part of the self energy is the same as that obtained 
in the study of neutrinos with standard model interactions in a medium in 49]. 

To leading order in perturbation theory the denominator in equation (|2.42p is SM^p. When the matter potential 
dominates (at high temperature in the standard model), 14a ^ SM"^ and SM'^p ^ Vaa oc G » ImSaa oc G^, 
thus in this regime 7 oc G <C 1. For example with active neutrinos with standard model interactions at high 
temperature, it was argued in ref. [i^ that Vaa GpkT^ /M^ whereas ImE^a ^ G^kT^ therefore at high temperature 
7 ~ ImSaa/Vaa ~ 5iu ^ 1 with g^, the standard model weak coupling. 

In the opposite limit, for SM'^ ^ Vaa oc G the vacuum mass difference dominates p ^ 1 and 7 <C 1 since SM"^ ^ 
G ^ G^. This analysis is similar to that in ref . [49| and precludes the possibility of "quantum zeno suppression" 0, [s^ 
at high temperature. 

The only region in which 7 may not be perturbatively small is near a resonance at which p ~ \ sin 20 1 and only for 
very small vacuum mixing angle so that 6M'^\ sin 20 1 cx G^. This situation requires a careful re-examination of the 
perturbative expansion, and in this case the propagator cannot be described as two separate Breit-Wigner resonances 
because the width of the resonances is of the same order of or larger than the separation between them. Such a 
possibility would require a complete re-assessment of the dynamics of the propagating modes in the medium as a 
consequence of the breakdown of the Breit-Wigner (or narrow width) approximation. However, for very small vacuum 
mixing angle, indeed a distinct possibility for keV sterile neutrinos 4^, the MSW resonance is very narrow and in 
most of the parameter range 7 <C 1 and can be safely neglected. This is certainly the case at very high or very low 
temperature regimes in which Vaa ^ SM^ or Vaa ^ respectively. 

In summary, it follows from this discussion that 7(fc) ^ 1, with the possible exception near an MSW resonance for 
extremely small vacuum mixing angle [48j. and such a case must be studied in detail non-perturbativcly. 

Hence, neglecting perturbatively small corrections, the Green's function in the flavor basis can be written as 



G(fc; t) = U{0^{k))G^{k; t)U'\e^{k)) (2.43) 

with 

/ g-i:i^t sin(Oi(fc)t) Q \ 

"^('^^ r^, 1 (2.44) 

\ " ^ n2(fe) / 

This Green's function and the expression for the damping rates ri^2 in eqn. (|2.36l2.37p lead to the following physical 
interpretation. The fields that diagonalize the Green's function on the mass shell, namely (f)i 2 are associated with 
the quasiparticle modes in the medium and describe the propagating excitations in the medium. From eqn. (|2.43|) 
these are related to the flavor fields (l)a,s by the unitary transformation 

00 = cos 01 -I- sin6',„(/)2 ; = cos6l„j02 - sinfl^^i . (2-45) 

When the matter potential Vaa ^ SAP, namely when 6m ~ 7r/2 it follows that (j)a ~ <p2 and the damping rate of the 
active species is r2 ~ Taa while 4>s ~ 4>i and the damping rate of the "sterile" species is Fi ~ Taa cos^ 9m Taa, 
where 

F.. 'JI^^^ (2.46) 
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is the ultrarelativistic limit of the damping rate of the active species in absence of mixing. In the opposite hmit, when 
the medium mixing angle is small 6^ ~ 0, corresponding to the near-vacuum case, ipa 4>i and the active species 
has a damping rate Fi ^ Taa while 0s ^ (j)2 with r2 ^ Taa sin^ 6,n ^ ^aa- In both limits the sterile species is weakly 
coupled to the plasma, active and sterile species become equally coupled near an MSW resonance for 6'm ~ 7''/4. 

We emphasize that the relation (|2.45p is not a relation between wave functions, but between the fields associated 
with the flavor eigenstates (active-sterile) and those associated with the propagating (quasiparticle) excitations in the 
medium (see section (fV| V 

III. QUANTUM KINETICS: 

The distribution functions for the active (a) and sterile (s) species are defined in terms of the diagonal entries of 
the mass matrix in the flavor representation, namely 



1 



{Mk; t)M-k, t)) + w^ik){Mk; t)M-k, t)) 



a 



a, s 



(3.1) 



where 



I^J(fc) = fc2 + Ml 



(3.2) 



The equal time expectation values of Heisenberg field operators are in the initial density matrix, and as shown in 
references m, [sj, Hs] they are the same as the equal time expectation value of the center of mass Wigner variables VP, 
where the expectation value is now in terms of the initial density matrix for the system and the distribution function 
of the noise which is determined by the thermal bath^S, >54„ Therefore the distribution functions for the active 
and sterile species are given by 



iV„(fc;i) 



1 



2Wo,ik) 



{i>a{k;t)^o.{~k,t))+Wl{k){^^{k;t)-i>^{-k,t)) 



(3.3) 



and the averages are taken over the initial density matrix of the system and the noise probability distribution. 
This expression combined with ean. (|2.28p makes manifest that the full time evolution of the distribution function is 
completely determined by the propagator Gap{k, t) obtained from the solution of the effective equations of motion in 
the medium [isj. 

It proves convenient to introduce a matrix of distribution functions in terms of a parameter VL as follows 



^I3{k,t; Vl) 



1 

2£l 



(*„(fc; i)*/3(-fc, t)) + r!2(*„(fc; t)*^(-fc, t)) 



1 



(3.4) 



from which we extract the active and sterile distribution functions from the diagonal elements, namely 

Na{k;t) = l^aAk.t;Wa{k)) ■ iV,(fc; t) - ]N,,,(A:, i; iy,(A:)) (3.5) 

and the off-diagonal elements determine off-diagonal correlation functions of the flelds and their canonical momenta 
in the flavor basis. 

We consider first the initial density matrix for the system p$(0) to be diagonal in the flavor basis with free field 
correlations 



1 



[l + 2N^{k;Q)]5a,f3 
[l + 2N^{k;0)] 6^p 



(3.6) 

(3.7) 
(3.8) 



with Na(k; 0) being the initial distribution functions for the active and sterile species. Different initial conditions will 
be studied below. 
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Following the steps described in appendix (|^ it is convenient to write ]N(fc,t;ri) = ¥(^^\k,t;D,) + ¥l^^\k,t;n) 
where IN*^^^ depends on the initial conditions but not on the noise f and IN'^^^ depends on the noise f but not on the 
initial conditions. We find 



+ lR,(2)(0^)g- 





cos2(0„O 








2Wan 






sin^(6'm) 








2Wan 





2 + ^a(O) 



+ sin2(6'„) 



21 r 



NaiO) 



1 

2 



2nw, 



2+^.(0) 



2nWa 



2nw, 



2nWs 

(Wa - Ws) ( ^1^2 



4f] \WaW, 



1 

2 

- 1 



29\ 



2^1 



(3.9) 



We have suppressed the dependence on k to simplify the notation. The contribution from the noise term can be 
written as 



]N(«)(fc,t;n) = ^ / ^[/(0„)|/i™(c^,i)/C™M/i:„(c^,t) + f^Vm(^,i)/C,„H/;;(a>,t)l[/-i(0™) 



where 



hm{oJ,t) = / e 
/o 



G^{k;t') ; U{uj,t)^ I e-'^' G^{k-t') 





and 



(3.10) 



(3.11) 



lC.m{uj) = ImEaa(fc; w) [1 + 2n(cj 

After lengthy but straightforward algebra we find 
]N(«Hfc,t;f7) = 



cos 



cos(6'm) sin(6'm) 



cos(6'„) sin(6'„) sin (6*™) 



"o2 + r!f 




\+n{n^{k)) 


(l-e-riW*) 


2r2if7 








\ + n{Sl2{k)) 


(l-e-^^W*) 


2Vt2VL 





(3.12) 



(3.13) 



where we have neglected terms proportional to 7. 

Approximations: In arriving at the expressions (|3.9p . p.l3p . we have made the following approximations: 

• (a) We have taken Z^jf' = 1 thus neglecting terms which are perturbatively small, of 0{G^). 

• (b) We have assumed Ti/^li <^ 1, which is warranted in perturbation theory and neglected terms proportional 
to this ratio. 

• (c) As discussed above, consistently with perturbation theory we have assumed 7(fc) ^ 1 and neglected terms 
proportional to it. This corresponds to the interaction rate much smaller than the oscillation frequencies and 
relies on the consistency of the perturbative expansion. 

• (d) In oscillatory terms we have taken a time average over the rapid time scales l/f2i^2 replacing sin^(f2i^20 = 
cos2(f]i_2t) ^ 1/2 ; sin(17i,2t) = cos(l]i,2i) ^0. 

Ultrarelativistic limit: The above expressions simplify considerably in the ultrarelativistic limit in which 

n ~ Wa{k) ~ w,{k) - ni{k) - n2{k) ~ fc, (3.14) 
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and in this limit it follows that 



Fl = Taa COS^ 9rn ] '^2 = Taa Sin^ 



(3.15) 



-rot 



and Taa is the ultrarelativistic limit of the width of the active species in the absence of mixing given by eqn. (j2.46p . 
In this limit we obtain the following simple expression for the time evolution of the occupation number matrix in the 
flavor basis (suppressing the k dependence for simplicity) 



m[t) = IR,(i)(0™) 7i(r!i)+(7Va(O)cos2(0„) + Ar,(O)sin2(0„)-n(r!i))e^i^i 

+ IR,('H^m) n{n^)+[Na{Q)sm^{e^) + NMcos\e,^)~n{n 
+ i R(3)(0^)g-i(ri+r.)t cos pi - r!2)i] {Na{Q) - N,{Q)) . 
It is straightforward to verify that 

The active and sterile populations are given by the diagonal elements of p.l6p . namely 
Nait) = cos^iOm) n(r!i)+(7Va(O)cos2(0„) + 7V,(O)sin2(0„)-n(r!i) 

+ sin2(0„) n(r!2)+ (iV,(O)sin2(0„) + 7V,(O)cos2(0™) -n(r!2)) e 
+ i sin2(20„)e-5(ri+r2)tcos[(fii-f72)t](A^a(O)-7V.(O) 



Fit 



(3.16) 



(3.17) 



(3.18) 



N,[t) = sin2(0„) 



+ (^Na{0) cos2(e„) + 7V,(0) sin2(^^„) - n(r!i) 



-Fit 



n(r!2) + (7V,(0) sin2((?„) + Ar,(0) cos''{9,n) - n(r!2)) 



I sin2(20„) 6-^(^1+^^)* cos [(1^1 - n2)t] (NaiO) - iV3(0)) . 



(3.19) 



The oscillatory term which results from the interference of the propagating modes 1, 2 damps out with a damping 
factor 



1 r 

2(ri + r2) 



aa 



(3.20) 



which determines the decoherence time scale Tdec = "^l^aa- These expressions are one of the main results of this 
article. 

Initial density matrix diagonal in the 1 — 2 basis: The above results were obtained assuming that the initial 
density matrix is diagonal in the flavor basis, if instead, it is diagonal in the basis of the propagating modes in the 
medium, namely the 1 — 2 basis, it is straightforward to find the result 



K{t) - u{e.^)< 



1 




i(f^i) + (a^i(0) -n(f^i) 




1 



i{n2) + [n2{q) - n{n2) 



(3.21) 
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In particular, the active and sterile distribution functions become 
Nait) = cos2 9,n + (^i(O) - n(r!i)) e-^^ 

N,{t) = cos^e^ n(r!2)+(A^2(0)-n(r!2))e 



-rat 



sm 



1(^(2)+ (iV2(0)-7l(r!2) 

(iVi(O) -n(17i) 



-Fit 



(3.22) 
(3.23) 



The results summarized by eqns. 



show that the distribution functions for the propagating modes in the 



medium, namely the 1, 2 quasiparticles, reach equilibrium with the damping factor ri,2 which is twice the damping 
rate of the quasiparticle modes (see eqn. (j2.44p V The interference term is present only when the initial density matrix 
is off diagonal in the (1,2) basis of propagating modes in the medium. 

If the initial density matrix is off-diagonal in the (1,2) basis, these off diagonal components damp-out within the 
decoherence time scale Tdec-, while the diagonal elements attain the values of the equilibrium distributions on the time 
scales l/ri,l/r2. 



IV. THE QUANTUM MASTER EQUATION 

The quantum master equation is the equation of motion of the reduced density matrix of the system fields in the 
interaction picture after integrating out the bath degrees of freedom. The first step is to define the interaction picture, 
for which a precise separation between the free and interaction parts in the Hamiltonian is needed [56]. In order to 
carry out the perturbative expansion in terms of the eigenstates in the medium, we include the lowest order forward 
scattering correction, namely the index of refraction into the un-perturbed Hamiltonian. This is achieved by writing 
the term 

4'lx' = ^l{x')+<l^l5x' (4.1) 

where 

Sx^^x^-ix') ; {Sx^) = o (4.2) 

and the average is performed in the bath density matrix /9_b(0) = e^^^"^^'^^ In this manner the quadratic part of 
the Lagrangian density for the active and sterile fields is 

£o[(f)] = i {^^^^a^* - 4>'^(M2 + V)$} (4.3) 

where V is the matter potential given by eqn. (|2.2ip . The unperturbed Hamiltonian for the system fields in the 
medium is diagonalized by the unitary transformation (|2.4p but with the unitary matrix U{Om) with 9m being the 
mixing angle in the medium given by equations (j2.33l2.34p and ipi , tp2 are now the fields associated with the eigenstates 
of the Hamiltonian in the medium including the index of refraction correction from the matter potential to 0{G) 
{0{Gf) in the case of neutrinos with standard model interactions). Introducing creation and annihilation operators 
for the fields 931,2 with usual canonical commutation relations, the unperturbed Hamiltonian /or the propagating modes 
in the medium including the index of refraction is 

HsM = E E ['^l{k)a^{kU{k)] (4.4) 

where Ui (fc) are the propagating frequencies in the medium given in equation (|2.38l2.39p . The interaction Hamiltonian 
is 

Hi = gJ d'x [4>l6x^ + 4>aOa] (4.5) 



where 



).^i 4- sin(e'™)(^2 • (4.6) 



This formulation represents a re-arrangement of the perturbative expansion in terms of the fields that create 
and annihilate the propagating modes in the medium. The remaining steps are available in the quantum optics 
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liter ature[5a|. Denoting the Hamiltonian for the bath degrees of freedom i?o[x, W^] = Hb the total Hamiltonian is 
H — Hg + Hb + Hj = Hq + Hj. The density matrix in the interaction picture is 

Mt) = e^"'p{t)e-''^°' (4.7) 
where p{t) is given by eqn. ()2.13|) and it obeys the equation of motion 

dpi{t) 



dt 



-i[Hi{t)Mt)] (4.8) 



with Hj{t) = e*^o*_ff/e ^^o* is the interaction Hamiltonian in the interaction picture of Hq. Iteration of this equation 
up to second order in the interaction yields [s^ 

^ = -^ [Hj{t),MO)] - £ dt' [Hj{t), [Hj{t'),Mt')]] + ' ' ' (4-9) 

The reduced density matrix for the system is obtained from the total density matrix by tracing over the bath 
degrees of freedom which are assumed to remain in equilibrium [Hg']. At this stage, several standard approximations 
are invoked [soj: 

• i): factorization: the total density matrix is assumed to factorize 

Mt) = PsAt) ® PsiO) (4.10) 

where it is assumed that the bath remains in equilibrium, this approximation is consistent with obtaining the 
effective action by tracing over the bath degrees of freedom with an equilibrium thermal density matrix. The 
correlation functions of the bath degrees of freedom are not modified by the coupling to the system. 

• ii): Markovian approximation: the memory of the evolution is neglected and in the double commutator in 
()4.9|) pi{t') is replaced by pi{t) and taken out of the integral. 

Taking the trace over the bath degrees of freedom yields the quantum master equation for the reduced density matrix, 

dps,i{i) 



dt 



dt'TvBi [Hiit), [Hi{t'), %{t)]] ] + ■■■ (4.11) 



where the first term has vanished by dint of the fact that the matter potential was absorbed into the unperturbed 
Hamiltonian, namely TrBPB{^)5'X^ — 0. This is an important aspect of the interaction picture in the basis of the 
propagating states in the medium. Up to second order we will only consider the interaction term 

Hiit) = J2 [cos0„</J,_j(<) +sin0„^2_^(t)] 0^j:{t) (4.12) 

k 

where we have written the interaction Hamiltonian in terms of spatial Fourier transforms and the fields are in the 
interaction picture of Hq. We neglect non-linearities from the second order contributions of the term (j>'^Sx'^, the non- 
linearities associated with the neutrino background are included in the forward scattering corrections accounted for 
in the matter potential. The quartic non-linearities are associated with active "neutrino-neutrino" elastic scattering 
and are not relevant for the production of the sterile species. 

The next steps are: i) writing out explicitly the nested commutator in (|4.1ip yielding four different terms, ii) taking 
the trace over the bath degrees of freedom yielding the correlation functions of the bath operators TrBO{t)0{t') 
(and t <-> t') and, iii) carrying out the integrals in the variable t' . While straightforward these steps are lengthy and 
technical and are relegated to appendix ([B|. Two further approximations are invoked [56|. 

• iii): the "rotating wave approximation": terms that feature rapidly varying phases of the form 
g±2iwi.2t. g±j(cji-hi^2)t g^j-g averaged out in time leading to their cancellation. This approximation also has a 
counterpart in the effective action approach in the averaging of rapidly varying terms, see the discussion after 
equation (|3.13p . 



iv): the Wigner Weisskopf approximation: time integrals of the form 



g-»(^-0)r^^ W -iV 







1 



- n 



ttS{uj - n) (4.13) 



where V stands for the principal part. The Markovian approximation (ii) when combined with the Wigner- 
Weisskopf approximation is equivalent to approximating the propagators by their narrow width Breit- Wigner 
form in the effective action. 
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All of these approximations i)- iv) detailed above are standard in the derivation of quantum master equations in 
the hterature|5a|. 

The quantum master equation is obtained in appendix ([BJ , it features diagonal and off-diagonal terms in the 1 — 2 
basis and is of the Lindblad form[56] which ensures that the trace of the reduced density matrix is a constant of 
motion as it must be, because it is consistently derived from the full Liouville evolution (|2.13[) . We now focus on the 
ultrarelativistic case uji{k) ^ uj2{k) ~ k which leads to substantial simplifications and is the relevant case for sterile 
neutrinos in the early Universe, we also neglect the second order corrections to the propagation frequencies. With 
these simplifications we obtain, 



dps,, 
dt 



where 



EE 



T,{k) 



1 + n{ujj{k))\ ps,ia^j{k)aj{k) + a^^{k)aj{k)ps,i - 2aj{k)ps,ia}j{k) 



+ n{ujj{k)) ps,taj{k)aUk) + aj{k)aUk)ps,t - 2a](k)ps,iaj{k) 



E 



m 



\ + n(uji{k))][ 4 (fc; t)ai (fc; t)ps,t + ps,ia[ (k\ t)a2 {k;t) - a2 {k; t)ps,ta\ {k; t) 



ai{k]t)ps,ia\{k;t)] + n{uji{k))[a2{k]t)a\{k]t)ps,i + ps,iai{k;t)a\{k;t) - a\{k;t)ps,iai{k]t) 



a\{k;t)ps,ia2{k;t) 



+ n{uj2{k))j (^a\{k] t)a2{k; t)ps,i + ps,ia\{k\ t)ai{k; t) ~ ai(fc; t)ps,ia\{k; t) 



a2{k;t)ps,ia\{k;t)) + n{uj2{k))[ai{k]t)a\{k;t)ps,i + ps,ia2{k;t)a\{k;t) - a\{k;t)ps,ia2{k;t) 



al{k;t)ps,iai{k;t) 



(4.14) 



ri(fc) =raa(fc)cOs2 



r(fc).-sin2.™r..(.) ; TUk)^ 



T2{k) = Taaik) sin^ 9m 
lmY,aaik,k) 



(4.15) 

and the interaction picture operators are given in eqn. (|B8[) . The expectation value of any system's operator A is 
given by 



{A){t) = TTp,^s{t)Ait) 



(4.16) 



where A(t) is the operator in the interaction picture of Hq, thus the time derivative of this expectation value contains 
two contributions 



-(A)(t) = TTp,,smit) + Trp,,5(t)A(t) , 



(4.17) 



The distribution functions for active and sterile species is defined as in equation (|3.3p with the averages defined as in 
(|4.16p . namely 



N^{k;t) =Trp,^sit) 



2Wa{k) 2 



4>a{k;t)(j)a{-k,t) 



(4.18) 



where the fields are in the interaction picture of Hq. The active and sterile fields are related to the fields that create 
and annihilate the propagating modes in the medium as 



(paik) = cos9„,^i{k) + sinemf2{k) ; (j)s{k) ^ cos 9,nf2{k) ~ sin e^fiik) . 
In the interaction picture of Hq 

1 



tpj{k,t) = 



,/%J~(k) 



,(fc) e-^^^C^)* +at(-fc) e''^' 



(k)t 



(4.19) 



(4.20) 
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where LOj{k) are the propagation frequencies in the medium up to leading order in G, given by equations (|2.38l2.39p . 
Introducing this expansion into the expression (|4.18p we encounter the ratio of the propagating frequencies in the 
medium ojj and the bare frequencies Wa- Just as we did in the previous section, we focus on the relevant case of 
ultrarelativistic species and approximate as in equation (|3.14p ujj{k) Wa{k) ~ k, in which case we find the relation 
between the creation- annihilation operators for the flavor fields and those of the 1, 2 fields to befisj 



aa{k, t) — cos 9mCii (fc, t) + siTi9mCi2ik, t) ; as(fc, t) — cos 0,„a2(fc, t) — sin 9mai{k, t) 
leading to the simpler expressions for the active and sterile distributions, 



Naik;t)^TTp,,sit) 



cos^ 9ma[{k, t)ai{k, t) + sin^ 9ma\{k, t)a2{k, t) 
1 



sin 2(?m(ai(fc, t)a2{k, t) + a\{k, t)ai{k, t)) 



(4.21) 



(4.22) 



N,{k-t)^TYp^,s{t) 



cos^ 9ma2{k, t)a2{k, t) + sin^ 0ma^(fc, t)ai{k, t) 



sin 29m{a\{k, t)a2{k, t) + al{k, t)ai (fc, <)) 



(4.23) 



In the interaction picture of Hq the products a-{k^t)aj{k^t) are time independent and a\{k,t)a2{k^t) — 
a| (fc, 0)a2(fc, 0) e«('^i('=)-'^2(fc))*. it is convenient to introduce the distribution functions and off-diagonal correlators 

nn(fc,<) =Trp,,s(t)at(fc,t)ai(fc,t) , n22{k,t) =i:Yp,,s{t)a\{k,t)a2{k,t) (4.24) 
ni2{k,t)=TTp,,s{t)a\{k,t)a2{k,t) , n2i{k,t) = i:Yp,^s{t)al(k,t)ai{k,t) ^ n\^{k,t) . (4.25) 

In terms of these, the distribution functions for the active and sterile species in the ultrarelativistic limit becomes 

Na{k]t) = cos2 6l,„nii(fc;t) -I- sin2 6l„n22(fc;t) + i sin 26l,„ (7112 (fc;t) 4-7121 (fc;t)) (4.26) 
Ns{k]t) = sin^ Till (fc; t) -f-cos^ 6l„ 7^22 (fc;t) - i sin 26l,„ (7112 (fc;t) 7121 (fc;t)) . (4.27) 

From eqn. (|4.17p we obtain the following kinetic equations for (fc; t) 

(4.28) 
(4.29) 
(4.30) 
(4.31) 



nil = 


-Fi 


^11 - UeqS 


r 

~ 2 


ni2 - 


1- 121 




-r2 


n22 - neqa 


r 

" 2 


ni2 - 


^21 



ni2 



"21 = 



— iAuj — 


r 1 

^ aa 

~Y. 


ni2 - 


r - 

" 2 . 


+ iAuj - 


r 1 

^ aa 

~Y. 


7121 - 


f - 
' 2 . 



(nil - 1eg,l) -I- (?122 - ?le,,2) 
(nil - legj) -f- (n22 - ne,^2) 



where n„ 



n{ujj{k)) are the equilibrium distribution functions for the corresponding propagating modes, and Aco = 



(a;2(fc) — ijJi(fc))- As we have argued above, in perturbation theory Taaik)/ ALij{k) <C 1, which is the same statement 
as the approximation 7 <C 1 as discussed for the effective action, and in this case the off diagonal contributions to the 
kinetic equations yield perturbative corrections to the distribution functions and correlators. To leading order in this 
ratio we find the distribution functions, 



nii{t) = neq,i + {nii[G) - rieg.i j e 
n22(i) = neg^2 + ("22(0) - neg,2 



Tit 



-Fit _ ^ 6 

2Aw 

-Fat 



2Aw 



i ni2(0) 
i ni2(0) 



+ c.c. 



+ c.c. 



(4.32) 
(4.33) 
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and ofF-diagonal correlators 



ni2(0 = e-'^^' e-—' {ni2{0) + i 



2Aw 



(nii(O) - neq,l) 



lAut g^(r2-ri)t _ j^' 



+ ("22(0) - neq,2) 



where 



r 1 ImSaJfc,fc) 1 . ~ 



2Aw 2 2fc(w2(fc) - wi(fc) 2 

with 7 defined in eqn. (|2.42[) and we have suppressed the momenta index for notational convenience. 



(4.34) 



(4.35) 



A. Comparing the effective action and quantum master equation 

We can now estabhsh the equivalence between the time evolution of the distribution functions obtained from the 
effective action and the quantum master equation, however in order to compare the results we must first determine 
the initial conditions in equations (j4. 32114. 34|) . The initial values 'T-ij(O) must be determined from the initial condition 
and depend on the initial density matrix. Two important cases stand out: i) an initial density matrix diagonal in the 
flavor basis or ii) diagonal in the 1 — 2 basis of propagating eigenstates in the medium. 

Initial density matrix diagonal in the flavor basis: the initial expectation values are obtained by inverting 
the relation between tpi^2 and 4'a,s- We obtain 

nii(O) - (al(fc)ai(fc))(0) = cos^ 0„,Ar,(O) + sin^ 0„iV,(O) (4.36) 
^22(0) = (aJ,(fc)a2(fc))(0) = cos2 0„,7V,(O) + sin^ 0„,iV„(O) (4.37) 

^12(0) - (al(fc)a2(fc))(0) = ^sm2e^{Na{0) - NM) (4-38) 

It is straightforward to establish the equivalence between the results obtained from the effective action and those 
obtained above from the quantum master equation as follows: i) neglect the second order frequency shifts (^^1,2 ~ ^1,2) 
and the perturbatively small corrections of order 7, ii) insert the initial conditions (|4. 36114.3^1) in the solutions (|4.321 
I4.34P , finally using the relations (|4.26l4.27p for the active and sterile distribution functions we find precisely the results 
given by equations (|3.18I3.19[) obtained via the non-equilibrium effective action. 

Initial density matrix diagonal in the 1 — 2 basis: in this case 

{a\{k)a,{k)){0) - 7Vi(0) 
(4(fc)a2(fc))(0) = 7V2(0) 

(al(fc)a2(fc))(0) = (4.39) 

with these initial conditions it is straightforward to obtain the result (|3.22l3.23p . 

The fundamental advantage in the method of the effective action is that it highlights that the main ingredient is 
the full propagator in the medium and the emerging time scales for the time evolution of distribution functions and 
coherences are completely determined by the quasiparticle dispersion relations and damping rates. 



B. Quantum kinetic equations: summary 



Having confirmed the validity of the kinetic equations via two independent but complementary methods, we now 
summarize the quantum kinetic equations in a form amenable to numerical study. For this purpose it is convenient 
to define the hermitian combinations 

riRit) = ni2{t) + n2i{t) ; ni{t) = i{ni2{t) ~ n2i{t)) (4.40) 

in terms of which the quantum kinetic equations for the distribution functions and coherences become (suppressing 
the momentum label) 
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Till = -TaaCOS 



"-22 



n-11 - nen.i 



1122 - neq.2 



■ sin 29m riR 



■ sin 29ra riR 



r r r 

flu = -{uj2 - uji) ni ^"^R |^sin26l,„ [nn ~ neq,i) + {n22 - neq,2) 

r 

ni = [uj2-^i)nR ^'^-f' 

with the active and sterile distribution functions related to the quantities above as follows 

Na{k] t) = cos^ 9mnii{k; t) + sin^ 0m"-22(fc; t) + \ sin29mnR{k; t) 



Ns{k]t) = sin 6',„nii(fc;i) +C0S 9„^n22{k;t) - - sm29mnR{k;t) 



(4.41) 
(4.42) 
(4.43) 
(4.44) 

(4.45) 
(4.46) 



In the perturbative limit when Taa sin 29 m/^uj <C 1 which as argued above is the correct limit in all but for a 
possible small region near an MSW resonance^], the set of kinetic equations simplify to 



nil 

?T'22 



-r„„ cos 9„ 



-1 aa Sm 



nil - nea.l 



n22 - ?^eq,2 



ur = - CJ2 - ^^1 ni 



ni 



-nR 



UJ2-L0l\ Ur 



-ni . 



In this case the active and sterile populations are given by (suppressing the momentum variable) 



(4.47) 
(4.48) 
(4.49) 
(4.50) 



Na{t) = COS2(0„O 



eq{uJl) + (nil(O) - UeqiuJl) 
3q{i^2) + (n-22(0) - neq{u)2) 



+ siv?{9m) 
+ sin(20m) e^~5^* cos [(wi - 072)^] ^^12(0) 



-T2t 



(4.51) 



N,{t) = s\n^{9m) 



COS^(0„ 



neqil^l) + (?^ll(0) - neg(t^l)^ 6 



7(^2) + (^?^22(0) - neq{u}2) 



-Tit 



-Fat 



where 



sin(26'm) e 2° * cos [{uji - U2)t] ni2(0) , 



ri(fc) =raa(fc)cos2 0,„ ; r2(fc) =r,,(A:)sin2 



(4.52) 



(4.53) 



and assumed that ni2(0) is real as is the case when the initial density matrix is diagonal both in the flavor or 1,2 
basis, 



(n\ — } 5 sin 20m {Na{0) — ^s(O)) diagonal in flavor basis 
' I diagonal in 1, 2 basis 



(4.54) 
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It is clear that the evolution of the active and sterile distribution functions cannot, in general, be written as simple 
rate equations. 

From the expressions given above for the quantum kinetic equations it is straightforward to generalize to account for 
the fermionic nature of neutrinos: the equilibrium distribution functions are replaced by the Fermi-Dirac distributions, 
and Pauli blocking effects enter in the explicit calculation of the damping rates. 



TRANSITION PROBABILITIES AND COHERENCES 



A "transition probability" in a medium 



The concept of a transition probability as typically used in neutrino oscillations is not suitable in a medium when 
the description is not in terms of wave functions but density matrices. However, an equivalent concept can be provided 
as follows. Consider expanding the active and sterile fields in terms of creation and annihilation operators. In the 
ultrarelativistic limit the positive frequency components are obtained from the relation (|4.2ip and their ensemble 
averages in the reduced density matrix are given by 



The kinetic equations for (ai.2(fc))(t) are found to be 

d 



dt 



'Mk)m 



zc.i(fc)-^^Vai(fc))(t) 



{a,{k)){t) 



Mk,t))-- (ai(fc))(0. 



where F has been defined in eqn. (|4.15p . To leading order in T/Aoj the solutions of these kinetic equations are 

if 



(ai(fc))(0 = (ai(A:))(0) e 



2Atj 
if 
2Atj 



(«2(fc))(0) 

(ai(fc))(0) 



(5.1) 

(5.2) 
(5.3) 

(5.4) 
(5.5) 



The initial values (ai,2(fc))(0) determine the initial values ipa,s{k; 0), or alternatively, giving the initial values ipa,sik; 0) 
determines (ai^2(fc))(0). Consider the case in which the initial density matrix is such that 



{aa{k,0)) = ipa{k)^0 ; (a,(fc, 0)) = (^,(fc, 0) = 

the initial values of (ai.2(fc))(0) are obtained by inverting the relation (14. 6p from which we find 

1 



Vs{k,t) 



■sin26'„(l - n) 



¥'a(A:,0) , 



this result coincides with that found in ref . [48[ . We can interpret the "transition probability" as 

2 



Pa->s[k,t) 



tpsik,t) 



fa{k,0) 



1 9 

- sin^ 20„ 



'Fit 



+ e"^^* - 2 cos ((w2 - e"^(^i+^^)* 



(5.6) 



(5.7) 



(5.8) 



where we have neglected perturbative corrections of 0{j). This result coincides with that obtain in ref.fS] from 
the effective action, and confirms a similar result for neutrinos with standard model interactions 49] . We emphasize 
that this "transition probability" is not obtained from the time evolution of single particle wave functions, but from 
ensemble averages in the reduced density matrix: the initial density matrix features a non- vanishing expectation value 
of the active field but a vanishing expectation value of the sterile field, however, upon time evolution the density 
matrix develops an expectation value of the sterile field. The relation between the transition probability (|5.8p and 
the time evolution of the distribution functions and coherences is now explicit, the first two terms in (|5.8[) precisely 
reflect the time evolution of the distribution functions ^11,7122 with time scales l/ri_2 respectively, while the last, 
oscillatory term is the interference between the active and sterile components and is damped out on the decoherence 
time scale Tdec- This analysis thus confirms the results in ref. [4^. 
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B. Coherences 

The time evolution of the off-diagonal coherence {a\{k)a2{k)){t) is determined by the kinetic equation (|4.30|) . 
neglecting perturbatively small corrections of 0(7) 

(al(fc)a2(fc))(i) - (al(fc)a2(fc))(0) e^^'^'e"^* (5.9) 

where we have used the relations (|4.15l ) in the ultrarelativistic limit. Therefore, in perturbation theory, if the initial 
density matrix is off-diagonal in the 1 — 2 basis (propagating modes in the medium) the off-diagonal correlations 
are exponentially damped out on the coherence time scale Tdec = ^/^aaik). This coherence term and its hermitian 
conjugate are precisely the ones responsible for the oscillatory term in the transition probability (jS.Sp . An important 
consequence of the damping of the off-diagonal coherences is that in perturbation theory the equilibrium de nsity 
matrix is diagonal in the basis of the propagating modes in the medium. This result confirms the arguments in ref. [5ll | . 
As can be seen from the expression of the transition probability (|5.8p this is precisely the time scale for suppression of 
the oscillatory interference term. However, the transition probability is not suppressed on this coherence time scale, 
the first two terms in (|5.8p reflect the fact that the occupation numbers build up on time scales l/Fi; l/r2 respectively 
and the interference term is exponentially suppressed on the decoherence time scale Tdec = 2/(ri + r2). For small 
mixing angle in the medium 9m all of these time scales can be widely different. 

It is noteworthy to compare the transition probability (|5.7p with the distribution functions (|4.51l4.52p . The first 
two, non-oscillatory terms in (j5.7p describe the same time evolution as the distribution functions 71-11,7122 of the 
propagating modes in the medium, while the last, oscillatory term describes the interference between these. This 
confirms the results and arguments provided in ref. [ioj . 



VI. FROM THE QUANTUM MASTER EQUATION TO THE QKE FOR THE "POLARIZATION" 

VECTOR 

The results of the previous section allows us to establish a correspondence between the quantum master equation 
(I4.14p the quantum kinetic equations (|4. 41114. 4^ and the quantum kinetic equation for a polarization vector often 
used in the literature[57„^]. Following ref. [59,], let us define the "polarization vector" with the following components. 



Po(fc,t) = (al{k,t)aaik,t)+al{k,t)asik,t)) =Naik,t)+N,{k,t) 

P,(fc,t) = {alik,t)a,ikA) + al{k,t)aaik,t)) 

Py{k,t) = -i{alik,t)as{k,t)-al{k,t)aaik,t)) 

P,{k,t) = {ai{k,t)aa{k,t)-a\{k,t)as{k,t))=Na{k,t)-Ns{k,t) 



(6.1) 
(6.2) 
(6.3) 
(6.4) 



where the creation and annihilation operators for the active and sterile fields are related to those that create and 
annihilate the propagating modes in the medium 1,2 by eqn. (|4.2ip . and the angular brackets denote expectation 
values in the reduced density matrix ps.i which obeys the quantum master equation (|4.14p . In terms of the population 
and coherences n.y the elements of the polarization vector are given by 



^0 = nil + 

Px = - sin 26m { nil - 7122 ) + cos 29„inR 



P, 



-ni 



Pz = COs20m^71ii - 7122^ + Sin20m77fl 

where nnj are defined by equation (|4.40p . Using the quantum kinetic equations (|4. 41114. 44|) we find 



dP^ 
dt 



-P. 



- aa 

"2~ 



{nil - 7leq,i^ + {n22 - neq,2^ 



■ COS 26r 



(6.5) 
(6.6) 
(6.7) 
(6.8) 



(6.9) 



dPx 
dt 



- (UJ2 — Wi) COS 29mni 



'- P 



■ sin 29r, 



T'eq,2 . 



(6.10) 
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LUi)nii 



(6.11) 



r r 

= -{u>2- u}i)sm20jnni —P^ — 



"-22 



^69,2 



We now approximate 



( \ (^2 -t^i) , . ^ „ 



thus neglecting the last terms in eqns. (|6.9l6.10p . introducing the vector V with components 

V—{uj2~uii) ( sin26'„, 0, — cos26'„ 



we find the following equations of motion for the polarization vector 



dp -f-t ^ Taa I 7-, - 

— = V X P PxX 

dt 2 ' 



dp 



» 



This equation is exactly of the form 



dt 



DPt 



dPo 

dt ' 



(6.12) 



(6.13) 



(6.14) 



(6.15) 



(6.16) 



used in the literature [Sg, |42, l43|, [571, l58|, where D and Pt can be identified from eqn. (I6.15p . 

Therefore the quantum kinetic equation for the polarization vector (|6.15p is equivalent to the full set of quantum 
kinetic equations (14. 41114. 44|) or equivalently to equations (|4.28II4.3"T|) under the approximation (|6.13p . Furthermore 
since the quantum kinetic equations (|4. 41114. 44| have been proven to be equivalent to the time evolution obtained from 
the effective action, we conclude that the kinetic equation for the polarization vector (|6.16p is completely equivalent 
to the effective action and the quantum master equation under the approximations discussed above. This equivalence 
between the effective action, the kinetic equations obtained from quantum Master equation and the kinetic equations 
for the polarization vector makes explicit that the fundamental scales for decoherence and damping are determined by 
ri^2, which are twice the damping rates of the quasiparticle modes. These are completely determined by the complex 
poles of the propagator in the medium. Furthermore the formulation in terms of the effective action, or equivalently 
the quantum master equation (|4.14p provides more information: for example from both we can extract the transition 
probability Pa—>s in the medium from expectation values of the field operators (or creation/annihilation operators) in 
the reduced density matrix, leading unequivocally to the expression (|5.8p which indeed features the two relevant time 
scales. Furthermore it directly yields information on the off-diagonal coherences (|5.9p which fall off on the decoherence 
time scale r^ec — '^/^aa, thus elucidating that the reduced density matrix in equilibrium (the asymptotic long time 
limit) is diagonal in the 1-2 basis. While this information could be extracted from linear combinations of Px,Py it 
is hidden in the solution of the kinetic equation for the polarization, whereas it is exhibited clearly in the quantum 
kinetic equations (|4. 28114. ST]) in the regime in which perturbation theory is applicable [Faa sin 20m/ (0^2 — ^ 1 ■ In 
this regime, which as argued above is the most relevant, the set of quantum kinetic equations (j4. 47114. 50]) combined 
with the relations (|4. 26114. 27)1 yield a much simpler and numerically amenable description of the time evolution of the 
populations and coherences: the active and sterile distribution functions are given by equations (|4.51l4.52p and the off- 
diagonal coherence by eqn. (|5.9p . Therefore, while the kinetic equation for the polarization and the quantum kinetic 
equations (|4. 47114. 50)) are equivalent and both are fundamentally consequences of the effective action or equivalently 
the quantum master equation, the study of sterile neutrino production in the early Universe does not implement any 
of these equivalent quantum kinetic formulations but instead assume a phenomenological approximate description in 
terms of a simple rate equation |4^ [43j . which implies only one damping scale. Such a simple rate equation cannot 
describe accurately the time evolution of distribution functions and coherences which involve two different time scales 
(away from MSW resonances). In our view, part of the problem in this formulation is the time averaged transition 
probability introduced in ref.[4^ which inputs the usual quantum mechanical vacuum transition probability but 
damped by a simple exponential on the decoherence time scale, clearly in contradiction with the result (j5.8p obtained 
from the reduced quantum density matrix. Within the kinetic formulation for the time evolution of the polarization 
vector Pq, eqn. (|6.15p it is not possible to extract the notion of a transition probability because the components of 
polarization vector are expectation values of bilinear operators in the reduced density matrix. Instead, the concept of 
active-sterile transition probability can be established in a medium via expectation values of the field operators (or 
their creation/annihilation operators) in the reduced density matrix es discussed in section (|V Ap . 
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VII. CONCLUSIONS 



Our goal is to study the non-equilibrium quantum kinetics of production of active and sterile neutrinos in a medium. 
We make progress towards that goal by studying a model of an active and a sterile mesons coupled to a bath in thermal 
equilibrium via couplings that model charged and current interactions of neutrinos. The dynamical aspects of mixing, 
oscillations, decoherence and damping are fairly robust and the results of the study can be simply modified to account 
for Pauli blocking effects of fermions and the detailed form of the matter potential. As already discussed in ref [48] with 
simple modifications, such as the detailed form of Vaa including the CP-odd and even terms(50|, and the Fermi-Dirac 
distributions for the equilibrium ones, this model provides a remarkably faithful description of the non-equilibrium 
dynamics of neutrinos. 

We obtained the quantum kinetic equations for the active and sterile species via two independent but complemen- 
tary methods. The first method obtains the non-equilibrium effective action for the active and sterile species after 
integrating out the bath degrees of freedom. This description provides a non-perturbative Dyson-like resummation 
of the self-energy radiative corrections, and the dynamics of the distribution functions is completely determined by 
the solutions of a Langevin equation with a noise term that obeys a generalized fluctuation-dissipation relation. The 
important ingredient in this description is the full propagator. The poles of the propagator correspond to two quasi- 
particle modes whose frequencies obey the usual dispersion relations of neutrinos in a medium with the corrections 
from the index of refraction (forward scattering), with damping rates (widths) 

Pi = Taa C0S2 0„ ; P2 = Paa sin^ . (7.1) 

where Paa is the interaction rate of the active species in absence of mixing (in the ultrarelativistic limit) and 9m, 
is the mixing angle in the medium. These two damping scales, along with the quasiparticle frequencies completely 
determine the evolution of the distribution functions. This is one of the important aspects of the kinetic description 
in terms of the non-equilibrium effective action: the dispersion relations and damping rates of the quasiparticle modes 
corresponding to the poles of the full propagator completely determine the non- equilibrium evolution of the distribution 
functions and coherences. 

We also obtained the quantum master equation for the reduced density matrix for the "neutrino degrees of freedom" 
by integrating (tracing) over the bath degrees of freedom taken to be in thermal equilibrium. An important aspect 
of the derivation consists in including the matter potential, or index of refraction from forward scattering to lowest 
order in the interactions in the unperturbed Hamiltonian. This method provides a re-arrangement of the perturbative 
expansion that includes self-consistently the index of refraction corrections and builds in the correct propagation 
frequencies in the medium. In this manner the the reduced density matrix (in the interaction picture ) evolves in 
time only through second order processes. From the reduced density matrix we obtain the quantum kinetic equations 
for the distribution functions and coherences. These are exactly the same as those obtained from the non- equilibrium 
effective action. We also obtain the kinetic equation for coherences and introduce a generalization of the active- 
sterile transition probability by obtaining the time evolution of expectation values of the active and sterile fields in 
the reduced quantum density matrix. Within the realm of validity of the perturbative expansion the set of kinetic 
equations for the distribution functions and coherences are given by 



fill = -Paa cos' 



2 , 



ni2 



nil - neq,l 
P 



"•22 = -Faa siu^ 6*,, 



"22 - n-e?,2 



{uj2{k) - uJi{k)^ Y ^^"^ ' "21 = "I2 ) (7-2) 



where "egj" = n{u!j{k)) are the equilibrium distribution functions for the corresponding propagating modes, wi,2(fc) 
are the dispersion relations in the medium including the index of refraction, and the active and sterile distribution 
functions are given by 

Na{k;t) = cos^e,nnii{k]t)+sm'^e„^n22{k;t) + ^sm2em(ni2{k;t) + n2i{k-,t)^ (7.3) 
Ns{k;t) = sm^emnii{k;t)+cos^9„in22{k;t)~^sm29m(^ni2{k;t) + n2i{k;t)y (7.4) 

The set of equations (17. 2p provide a simple system of uncoupled rate equations amenable to numerical study, 
whose solution yields the active and sterile distribution functions via the relations (|7.3l7.4p . with straightforward 
modifications for fermions. 
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From the kinetic equations above, it is found that the coherences 

"12 = {a\a2) (7.5) 

which are off-diagonal (in the 1 — 2 basis of propagating modes in the medium) expectation values in the reduced 
quantum density matrix are exponentially suppressed on a decoherence time scale r^ec = ^/Taa indicating that the 
equilibrium reduced density matrix is diagonal in the 1 — 2 basis, confirming the arguments in ref. 51]. 

The generalization of the active-sterile transition probability in the medium via the expectation value of the active 
and sterile fields in the reduced quantum density matrix yields 



Pa^s = ^ sin^ 20„ 



e"^i*-fe"^^*-2cos 



(^(UJ2 - UJi)t 



(7.6) 



this result shows that the active-sterile transition probability depends on the two damping time scales of the quasi- 
particle modes in the medium which are also the time scales of kinetic evolution of the distribution functions, and 
confirms the results of refs.|4^. 

Finally, from the full set of quantum kinetic equations (|4. 47114. 50|) and the approximation (|6.13p we have obtained 
the set of quantum kinetic equations for the polarization vector, most often used in the literature, 

where the relation between the components of the polarization vector Pq , P and the distribution functions and coher- 
ences is explicitly given by eqns. (|6.m6.4p (or equivalently (|6.5ll6.8p ). and V is given by eqn. (|6.14p . Thus we have 
unambiguously established the direct relations between the effective action, quantum master equation, the full set 
of kinetic equations for population and coherences and the quantum kinetic equations in terms of the "polarization 
vector" most often used in the literature. These are all equivalent, but the effective action approach distinctly shows 
that the two independent fundamental damping scales are those associated with Fi 2, namely the damping rates of 
the quasiparticles in the medium, which are determined by the complex poles of the propagator. Furthermore in the 
regime of validity of perturbation theory, the set of kinetic equations (|7.2p obtained from the quantum master equation 
yield a simple and clear understanding of the different time scales for the active and sterile distribution functions and 
a remarkably concise description of active and sterile production when combined with the relations (|7.3l7.4p . These 
simpler set of rate equations are hidden in the kinetic equation (|7.7p . 

We have also argued that the simple phenomenological rate equation used in numerical studies of sterile neu- 
trino production in the early Universe is not an accurate description of the non-equilibrium evolution, and trace its 
shortcomings to the time integral of an overly simplified description of the transition probability in the medium. 

Our study focused on a scalar model that features many similarities to but also distinct differences with the theory 
of mixed neutrinos. The dispersion relations, medium dependence of the mixing angles, transition probabilities for 
ensemble averages, and dependence of the damping rates of the propagating modes on the active collision rate are 
robust features in common with the case of neutrinos. These similarities are strengthened by the fact that the kinetic 
equations obtained in this article are identical to those available in the literature in terms of the polarization vector, 
with the bonus that we provide a different interpretation that highlights the role of the non-equilibrium evolution in 
terms of the physical propagating modes. All of these similarities and the combination of results obtained in this 
study and those reported in[48, 49] lend support to the expectation that the results obtained in this study are relevant 
for the description of the kinetics of neutrinos. 

There are, however, differences with the neutrino case that must eventually be addressed for a more complete 
treatment and understanding: spinorial and chiral structures, although these are not directly accounted for either in 
the quantum mechanical description of neutrino oscillations nor in the phenomenological description of the kinetics, 
Fermionic nature of the neutrino field, which enters in the distribution function, however, the simplicity of the kinetic 
equations found in this article allow a simple replacement of the distribution functions by the Fermi-Dirac one, 
automatically including Pauli blocking, furthermore, the matter potential in the case of neutrinos features both a CP- 
odd term arising from the lepton and baryon asymmetry, and a CP-even term that depends solely on temperature, 
the overall sign of the matter potential is determined by these two contributions. For the case of sterile neutrinos 
with keV masses, an MSW resonance is only available when the CP-odd term dominates. Our study in this article is 
general, without specifying a particular form of the matter potential and addressed all possibilities with or without 
MSW resonances. The only specific aspect is that the matter potential is flavor diagonal and only features an entry 
in the active- active matrix element. 

While the model studied here is clearly a simplification of the case of neutrinos, the body of results and similarities 
established with the neutrino case suggest a reliable description of the quantum kinetics. A more detailed study of 
the impact of the differences on the non-equilibrium dynamics will be the subject of forthcoming work. 
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APPENDIX A: A SIMPLER CASE 



Consider for simplicity the case of one scalar field. The solution of the Langevin equation is given by 

')dt' , (Al) 



^^{t) = g{k; t) + g{k; t) H^, + f g{k; t') ^^{t - t' 

Jo 



where the dot stands for derivative with respect to time. In the Breit-Wigner approximation and setting = 1 



n{k) 

where is the position of the quasiparticle pole (dispersion relation) and its width is given by 

^i{n{k)) 



m 



n{k) 



The particle number is given by 
N{k,t) = 



2W{k) 



(^'(fc, t)*(~fc, t)} + VF2(fc)(*(fc, t)} 



(A2) 



(A3) 



(A4) 



where W{k) is the bare frequency. Taking the initial density matrix of the field ^' to be that corresponding to a 
free-field with arbitrary non-equilibrium initial distribution function N(k; 0) and carrying out both averages, over the 
initial density matrix for the field and of the quantum noise and using that the average of the latter vanishes, we find 



N{k;t)^Niik;t)+N2{k;t)~^ 



with 



Niik;t) 



l + 27V(fc;0) 



AW^{k) 



(ff(fc; t)f + 2W\k){g(k; t)f + W\k)g\k- 1) 



W^{k)\h{ujM+\f{uoM 



where 



h{uj,t) 



-iujt 



-iujt' 



g{k;t')dt' 



g{k;t')dt' 



(A5) 

(A6) 
(A7) 

(A8) 
(A9) 

(AlO) 



The terms Ni{k;t); N2{k;t) have very different origins: the term Ni{k;t) depends on the initial condition and origi- 
nates in the first two terms in Ij2.28p namely those independent of the noise, which survive upon taking the average 
over the noise. The term N2{k\ t) is independent of the initial conditions and is solely determined by the correlation 
function of the noise term and is a consequence of the fluctuation dissipation relation. Using the expression (jA2p we 
find 



Ni{k-t) 



1 



7V(fc;0) 



-T{k)t 



1 + sin2(Sl(/c)t) 



Q?{k)-W^{k) 
2W{k)VL{k) 



O 



f^2(fc) 



(All) 



25 



where the neglected terms of order r^(fc)/fl^(fc) ^ 1 are perturbatively smalL The oscillatory term in (jAlip averages 
out on a short time scale l/il{k) ^ l/r(fc) and we can replace (jAlip by its average over this short time scale yielding 



Ni{k;t) 



-+iV(fc;0) 



-r{k)t 



1 f n^k)-w^k) 

^ ^ 2 2W{k)n{k) 



o 



(A12) 



In perturbation theory ^'^{k)-W'^{k)/2W{k)^{k) < 1, can be neglected to leading order in perturbative quantities, 
thus we obtain 



Ni{k;t) 



N{k;0) 



-r(k)t 



(A13) 



Using the fact that — — tj) we can perform the integrals in N2{k;t) in the narrow width (Breit-Wigner) 

approximation by using eqn. (|A3p . with the result 



N,ik;t)^Z^ 
Replacing in perturbation theory 



w^ik) + n^{k) 



2W{k)nik) 



Zk « 1 



n{n{k)) 



wHk) + n^k) 

2W{k)n{k) 



we find 

N{k; t) = n{n{k)) + {N{k; 0) - n{n{k))) e'^^'^)* 
which is the solution of the usual kinetic equation 

dN{k;t) 



dt 



where 



-r(fc) {N{k-t)-N,g{k)) 



N,g{k) n{n{k)) 



(AM) 

(A15) 

(A16) 

(A17) 
(A18) 



It is important to highlight the series of approximations that led to this result: i) the narrow width (Breit-Wigner) 
approximation, ii) Zk ^ 1, iii) fl^{k) ~ W'^{k), iv) T{k)/il{k) <^ 1, these approximations are all warranted in 
perturbation theory. Clearly including perturbative corrections lead to perturbative departures of the usual kinetic 
equation and of the equilibrium distribution function. 



APPENDIX B: QUANTUM MASTER EQUATION 



Taking the trace over the bath variables with the factorized density matrix (j4.10p . the double commutator in 
equation (|4.1ip becomes 



^PsMMt')Mt) TTBPBio)oit')oit) 

^Mt)PsAt)Mt') TrBPB{0)O{t')O{t) 
-Mt')psAt)Mt) TrBPi3(0)O(t)O(<')} 



(Bl) 



We suppressed the momentum index to simplify notation but used the fact that translational invariance of the bath 
implies that the correlation functions are diagonal in momentum. The bath correlation functions were given in ref. j48j 
(see section 3-B in this reference) and we just summarize these results: 



TTBPB{0)O{t)O{t') 
TTBPB{0)O{t')O{t) 



duo ImSaa(fc;w)[l+n(w)] g-'''^^*-*') 

duo ImSaa(fc;w) n{uj) e^^'^^*-*') 



(B2) 
(B3) 
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where we used the property ImEQa(fc; w) = —lmY}aa{k; — a;)[43|. The self energy E is obtained from the discontinuity 
across the W — x hues in the diagram in fig. Q and is the same quantity that enters in the non-equihbrium effective 
action, and n{Lo) is the equihbrium distribution function. The active field 0a is related to the fields that create and 
annihilate the propagating modes in the medium Lpi 2 as in eq. (|4.19p . hence terms of the form 



1 



(B4) 



and all other terms in (jBf p are written accordingly. The next step requires writing these fields in terms of creation 
and annihilation operators in the interaction picture of Hq^ their expansion is shown in eqn. (I4.20|) . The resulting 
products of creation and annihilation operators all feature phases which are re-arranged to depend separately on the 
variable t and t — t', for example 



^2iujjt ~iujj{t — t ) 



(t~t') 



, etc. . 



(B5) 



The exponentials that depend on t — t' , such as e^^^^^* * •* are combined with the exponentials in (|B2IB3[) and the 
integral in t' in (jBip is written as an integral in r = t — t'. The Wigner-Weisskopf approximation for the resulting 
integral yields eqn. (|4.13p . After performing the time integral the terms of the form a^aj do not feature any phase, 
whereas terms of the form aiUj (and their hermitian conjugate) feature terms of the form e='='('^»+"j)*, all of these 
rapidly oscillating terms average out and are neglected in the "rotating wave approximation" [s^ , which is tantamount 
to time-averaging these rapidly varying terms. The remaining terms can be gathered together into two different type 
of contributions, diagonal and off-diagonal in the 1 — 2 indices. The diagonal contributions do not feature explicit 
time dependence while the off-diagonal one features an explicit time dependence of the form "'^2)* _ 

Diagonal: The diagonal contributions are 



dps,, 
dt 



Y.T.\~ Wk)a,{kl PS At)] - ^ (1 + n{uj,{k))) (ps.,a]{k)a,{k) + a]{k)a,{k)P^ 



2aj{k)ps,ia]{k) I + n{ujj{k)) ( ps,iaj{k)a]{k) + aj{k)a]{k)ps,t - 2a]{k)ps,iaj(k) 



(B6) 



where the second order frequency shifts Awj (k) and the widths Tj (k) are given in equations 

Off diagonal: The full expression for the off-diagonal contributions is lengthy and cumbersome and we just quote 
the result for the real part of the quantum master equation, neglecting the imaginary part which describes a second 
order shift to the oscillation frequencies of the off-diagonal coherences. 



dps,-, 
dt 



ri(fc) 



(^1 + n{LUi{k)fj {a\{k]t)ai{k;t)ps,i + ps,ia\{k\t)a2{k;t) - a2{k;t)ps,ia\{k]t) 



ai (k; t)ps,ia\{k, t)j + n(cji (fc)) (^a2{k; t)a\ (fc; t)ps,i + ps^idi {k; t)a\{k- 1) - a\{k- t)ps,iai (fc; t) 



a\{k;t)ps,ia2{k;t) 



r2(fc) 



(^l + n{u!2{k))j (^al{k; t)a2{k; t)ps,t + ps,ia2{k;t)ai{k; t) - ai(fc; t)ps,ia2{k;t) 



a2{k;t)ps,ia\{k,t)j + n{ui2{k)){ai{k;t)a\(k,t)ps,i + ps,ia2{k;t)a\{k;t) - a\{k;t)ps,ia2(k,t) 



al{k;t)ps,iai{k;t) 



(B7) 



where the interaction picture operators 



aj{k; t) = aj{k\ 0) e 



(B8) 



and 



2 ^J uji{k)u)2{k) 



(B9) 
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